Large-scale phosphoproteomics reveals activation of the MAPK/GADD45ß/P38 axis and cell cycle inhibition in response to BMP9 and BMP10 stimulation in endothelial cells.

BACKGROUND: BMP9 and BMP10 are two major regulators of vascular homeostasis. These two ligands bind with high affinity to the endothelial type I kinase receptor ALK1, together with a type II receptor, leading to the direct phosphorylation of the SMAD transcription factors. Apart from this canonical pathway, little is known. Interestingly, mutations in this signaling pathway have been identified in two rare cardiovascular diseases, hereditary hemorrhagic telangiectasia and pulmonary arterial hypertension. METHODS: To get an overview of the signaling pathways modulated by BMP9 and BMP10 stimulation in endothelial cells, we employed an unbiased phosphoproteomic-based strategy. Identified phosphosites were validated by western blot analysis and regulated targets by RT-qPCR. Cell cycle analysis was analyzed by flow cytometry. RESULTS: Large-scale phosphoproteomics revealed that BMP9 and BMP10 treatment induced a very similar phosphoproteomic profile. These BMPs activated a non-canonical transcriptional SMAD-dependent MAPK pathway (MEKK4/P38). We were able to validate this signaling pathway and demonstrated that this activation required the expression of the protein GADD45ß. In turn, activated P38 phosphorylated the heat shock protein HSP27 and the endocytosis protein Eps15 (EGF receptor pathway substrate), and regulated the expression of specific genes (E-selectin, hyaluronan synthase 2 and cyclooxygenase 2). This study also highlighted the modulation in phosphorylation of proteins involved in transcriptional regulation (phosphorylation of the endothelial transcription factor ERG) and cell cycle inhibition (CDK4/6 pathway). Accordingly, we found that BMP10 induced a G1 cell cycle arrest and inhibited the mRNA expression of E2F2, cyclinD1 and cyclinA1. CONCLUSIONS: Overall, our phosphoproteomic screen identified numerous proteins whose phosphorylation state is impacted by BMP9 and BMP10 treatment, paving the way for a better understanding of the molecular mechanisms regulated by BMP signaling in vascular diseases.

Large-Scale Chromosomal Changes Lead to Genome-Level Expression Alterations, Environmental Adaptation, and Speciation in the Gayal (Bos frontalis).

Determining the functional consequences of karyotypic changes is invariably challenging because evolution tends to obscure many of its own footprints, such as accumulated mutations, recombination events, and demographic perturbations. Here, we describe the assembly of a chromosome-level reference genome of the gayal (Bos frontalis) thereby revealing the structure, at base-pair-level resolution, of a telo/acrocentric-to-telo/acrocentric Robertsonian translocation (2;28) (T/A-to-T/A rob[2;28]). The absence of any reduction in the recombination rate or genetic introgression within the fusion region of gayal served to challenge the long-standing view of a role for fusion-induced meiotic dysfunction in speciation. The disproportionate increase noted in the distant interactions across pro-chr2 and pro-chr28, and the change in open-chromatin accessibility following rob(2;28), may, however, have led to the various gene expression irregularities observed in the gayal. Indeed, we found that many muscle-related genes, located synthetically on pro-chr2 and pro-chr28, exhibited significant changes in expression. This, combined with genome-scale structural variants and expression alterations in genes involved in myofibril composition, may have driven the rapid sarcomere adaptation of gayal to its rugged mountain habitat. Our findings not only suggest that large-scale chromosomal changes can lead to alterations in genome-level expression, thereby promoting both adaptation and speciation, but also illuminate novel avenues for studying the relationship between karyotype evolution and speciation.

Different cardiovascular and pulmonary phenotypes for single- and double-knock-out mice deficient in BMP9 and BMP10.

AIMS: BMP9 and BMP10 mutations were recently identified in patients with pulmonary arterial hypertension, but their specific roles in the pathogenesis of the disease are still unclear. We aimed to study the roles of BMP9 and BMP10 in cardiovascular homeostasis and pulmonary hypertension using transgenic mouse models deficient in Bmp9 and/or Bmp10. METHODS AND RESULTS: Single- and double-knockout mice for Bmp9 (constitutive) and/or Bmp10 (tamoxifen inducible) were generated. Single-knock-out (KO) mice developed no obvious age-dependent phenotype when compared with their wild-type littermates. However, combined deficiency in Bmp9 and Bmp10 led to vascular defects resulting in a decrease in peripheral vascular resistance and blood pressure and the progressive development of high-output heart failure and pulmonary hemosiderosis. RNAseq analysis of the lungs of the double-KO mice revealed differential expression of genes involved in inflammation and vascular homeostasis. We next challenged these mice to chronic hypoxia. After 3 weeks of hypoxic exposure, Bmp10-cKO mice showed an enlarged heart. However, although genetic deletion of Bmp9 in the single- and double-KO mice attenuated the muscularization of pulmonary arterioles induced by chronic hypoxia, we observed no differences in Bmp10-cKO mice. Consistent with these results, endothelin-1 levels were significantly reduced in Bmp9 deficient mice but not Bmp10-cKO mice. Furthermore, the effects of BMP9 on vasoconstriction were inhibited by bosentan, an endothelin receptor antagonist, in a chick chorioallantoic membrane assay. CONCLUSIONS: Our data show redundant roles for BMP9 and BMP10 in cardiovascular homeostasis under normoxic conditions (only combined deletion of both Bmp9 and Bmp10 was associated with severe defects) but highlight specific roles under chronic hypoxic conditions. We obtained evidence that BMP9 contributes to chronic hypoxia-induced pulmonary vascular remodelling, whereas BMP10 plays a role in hypoxia-induced cardiac remodelling in mice.

UV-induced reduction in Polycomb repression promotes epidermal pigmentation.

Ultraviolet (UV) radiation is a prime environmental stressor that our epidermis is exposed to on a daily basis. To avert UV-induced damage, epidermal stem cells (EpSCs) become pigmented via a process of heterotypic interaction between melanocytes and EpSCs; however, the molecular mechanisms of this interaction are not well understood. In this study, we show that the function of a key chromatin regulator, the Polycomb complex, was reduced upon UV exposure in human and mouse epidermis. Genetic ablation of key Polycomb subunits in murine EpSCs, mimicking depletion upon UV exposure, results in an increased number of epidermal melanocytes and subsequent epidermal pigmentation. Genome-wide transcriptional and chromatin studies show that Polycomb regulates the expression of UV-responsive genes and identifies type II collagen (COL2A1) as a critical secreted regulator of melanogenesis and epidermal pigmentation. Together, our findings show how UV exposure induces Polycomb-mediated changes in EpSCs to affect melanocyte behavior and promote epidermal pigmentation.

Genomes reveal selective sweeps in kiang and donkey for high-altitude adaptation.

Over the last several hundred years, donkeys have adapted to high-altitude conditions on the Tibetan Plateau. Interestingly, the kiang, a closely related equid species, also inhabits this region. Previous reports have demonstrated the importance of specific genes and adaptive introgression in divergent lineages for adaptation to hypoxic conditions on the Tibetan Plateau. Here, we assessed whether donkeys and kiangs adapted to the Tibetan Plateau via the same or different biological pathways and whether adaptive introgression has occurred. We assembled a de novo genome from a kiang individual and analyzed the genomes of five kiangs and 93 donkeys (including 24 from the Tibetan Plateau). Our analyses suggested the existence of a strong hard selective sweep at the EPAS1 locus in kiangs. In Tibetan donkeys, however, another gene, i.e., EGLN1, was likely involved in their adaptation to high altitude. In addition, admixture analysis found no evidence for interspecific gene flow between kiangs and Tibetan donkeys. Our findings indicate that despite the short evolutionary time scale since the arrival of donkeys on the Tibetan Plateau, as well as the existence of a closely related species already adapted to hypoxia, Tibetan donkeys did not acquire adaptation via admixture but instead evolved adaptations via a different biological pathway.

Polycomb complexes redundantly maintain epidermal stem cell identity during development.

Polycomb repressive complex 1 (PRC1) and PRC2 are critical epigenetic developmental regulators. PRC1 and PRC2 largely overlap in their genomic binding and cooperate to establish repressive chromatin domains demarcated by H2AK119ub and H3K27me3. However, the functional contribution of each complex to gene repression has been a subject of debate, and understanding of its physiological significance requires further studies. Here, using the developing murine epidermis as a paradigm, we uncovered a previously unappreciated functional redundancy between Polycomb complexes. Coablation of PRC1 and PRC2 in embryonic epidermal progenitors resulted in severe defects in epidermal stratification, a phenotype not observed in the single PRC1-null or PRC2-null epidermis. Molecular dissection indicated a loss of epidermal identity that was coupled to a strong derepression of nonlineage transcription factors, otherwise repressed by either PRC1 or PRC2 in the absence of its counterpart. Ectopic expression of subsets of PRC1/2-repressed nonepidermal transcription factors in wild-type epidermal stem cells was sufficient to suppress epidermal identity genes, highlighting the importance of functional redundancy between PRC1 and PRC2. Altogether, our studies show how PRC1 and PRC2 function as two independent counterparts, thereby providing a repressive safety net that protects and preserves lineage identity.

Transcriptome analysis of neural progenitor cells derived from Lowe syndrome induced pluripotent stem cells: identification of candidate genes for the neurodevelopmental and eye manifestations.

BACKGROUND: Lowe syndrome (LS) is caused by loss-of-function mutations in the X-linked gene OCRL, which codes for an inositol polyphosphate 5-phosphatase that plays a key role in endosome recycling, clathrin-coated pit formation, and actin polymerization. It is characterized by congenital cataracts, intellectual and developmental disability, and renal proximal tubular dysfunction. Patients are also at high risk for developing glaucoma and seizures. We recently developed induced pluripotent stem cell (iPSC) lines from three patients with LS who have hypomorphic variants affecting the 3' end of the gene, and their neurotypical brothers to serve as controls. METHODS: In this study, we used RNA sequencing (RNA-seq) to obtain transcriptome profiles in LS and control neural progenitor cells (NPCs). RESULTS: In a comparison of the patient and control NPCs (n = 3), we found 16 differentially expressed genes (DEGs) at the multiple test adjusted p value (padj) < 0.1, with nine at padj < 0.05. Using nominal p value < 0.05, 319 DEGs were detected. The relatively small number of DEGs could be due to the fact that OCRL is not a transcription factor per se, although it could have secondary effects on gene expression through several different mechanisms. Although the number of DEGs passing multiple test correction was small, those that were found are quite consistent with some of the known molecular effects of OCRL protein, and the clinical manifestations of LS. Furthermore, using gene set enrichment analysis (GSEA), we found that genes increased expression in the patient NPCs showed enrichments of several gene ontology (GO) terms (false discovery rate < 0.25): telencephalon development, pallium development, NPC proliferation, and cortex development, which are consistent with a condition characterized by intellectual disabilities and psychiatric manifestations. In addition, a significant enrichment among the nominal DEGs for genes implicated in autism spectrum disorder (ASD) was found (e.g., AFF2, DNER, DPP6, DPP10, RELN, CACNA1C), as well as several that are strong candidate genes for the development of eye problems found in LS, including glaucoma. The most notable example is EFEMP1, a well-known candidate gene for glaucoma and other eye pathologies. CONCLUSION: Overall, the RNA-seq findings present several candidate genes that could help explain the underlying basis for the neurodevelopmental and eye problems seen in boys with LS.

Potential dual expansion of domesticated donkeys revealed by worldwide analysis on mitochondrial sequences.

Molecular studies on donkey mitochondrial sequences have clearly defined two distinct maternal lineages involved in domestication. However, domestication histories of these two lineages remain enigmatic. We therefore compared several population characteristics between these two lineages based on global sampling, which included 171 sequences obtained in this study (including Middle Asian, East Asian, and African samples) plus 536 published sequences (including European, Asian, and African samples). The two lineages were clearly separated from each other based on whole mitochondrial genomes and partial non-coding displacement loop (D-loop) sequences, respectively. The Clade I lineage experienced an increase in population size more than 8 000 years ago and shows a complex haplotype network. In contrast, the population size of the Clade II lineage has remained relatively constant, with a simpler haplotype network. Although the distribution of the two lineages was almost equal across the Eurasian mainland, they still presented discernible but complex geographic bias in most parts of Africa, which are known as their domestication sites. Donkeys from sub-Saharan Africa tended to descend from the Clade I lineage, whereas the Clade II lineage was dominant along the East and North coasts of Africa. Furthermore, the migration routes inferred from diversity decay suggested different expansion across China between the two lineages. Altogether, these differences indicated non-simultaneous domestication of the two lineages, which was possibly influenced by the response of pastoralists to the desertification of the Sahara and by the social expansion and trade of ancient humans in Northeast Africa, respectively.

Convergent genomic signatures of high-altitude adaptation among domestic mammals.

Abundant and diverse domestic mammals living on the Tibetan Plateau provide useful materials for investigating adaptive evolution and genetic convergence. Here, we used 327 genomes from horses, sheep, goats, cattle, pigs and dogs living at both high and low altitudes, including 73 genomes generated for this study, to disentangle the genetic mechanisms underlying local adaptation of domestic mammals. Although molecular convergence is comparatively rare at the DNA sequence level, we found convergent signature of positive selection at the gene level, particularly the EPAS1 gene in these Tibetan domestic mammals. We also reported a potential function in response to hypoxia for the gene C10orf67, which underwent positive selection in three of the domestic mammals. Our data provide an insight into adaptive evolution of high-altitude domestic mammals, and should facilitate the search for additional novel genes involved in the hypoxia response pathway.

PHF19 promotes multiple myeloma tumorigenicity through PRC2 activation and broad H3K27me3 domain formation.

Dysregulation of polycomb repressive complex 2 (PRC2) promotes oncogenesis partly through its enzymatic function for inducing trimethylation of histone H3 lysine 27 (H3K27me3). However, it remains to be determined how PRC2 activity is regulated in normal and diseased settings. We here report a PRC2-associated cofactor, PHD finger protein 19 (PHF19; also known as polycomb-like 3), as a crucial mediator of tumorigenicity in multiple myeloma (MM). Overexpression and/or genomic amplification of PHF19 is found associated with malignant progression of MM and plasma cell leukemia, correlating to worse treatment outcomes. Using various MM models, we demonstrated a critical requirement of PHF19 for tumor growth in vitro and in vivo. Mechanistically, PHF19-mediated oncogenic effect relies on its PRC2-interacting and chromatin-binding functions. Chromatin immunoprecipitation followed by sequencing profiling showed a critical role for PHF19 in maintaining the H3K27me3 landscape. PHF19 depletion led to loss of broad H3K27me3 domains, possibly due to impaired H3K27me3 spreading from cytosine guanine dinucleotide islands, which is reminiscent to the reported effect of an "onco"-histone mutation, H3K27 to methionine (H3K27M). RNA-sequencing-based transcriptome profiling in MM lines also demonstrated a requirement of PHF19 for optimal silencing of PRC2 targets, which include cell cycle inhibitors and interferon-JAK-STAT signaling genes critically involved in tumor suppression. Correlation studies using patient sample data sets further support a clinical relevance of the PHF19-regulated pathways. Lastly, we show that MM cells are generally sensitive to PRC2 inhibitors. Collectively, this study demonstrates that PHF19 promotes MM tumorigenesis through enhancing H3K27me3 deposition and PRC2's gene-regulatory functions, lending support for PRC2 blockade as a means for MM therapeutics.

Large-scale phylogenetic analyses provide insights into unrecognized diversity and historical biogeography of Asian leaf-litter frogs, genus Leptolalax (Anura: Megophryidae).

Southeast Asia and southern China (SEA-SC) harbor a highly diverse and endemic flora and fauna that is under increasing threat. An understanding of the biogeographical history and drivers of this diversity is lacking, especially in some of the most diverse and threatened groups. The Asian leaf-litter frog genus Leptolalax Dubois 1980 is a forest-dependent genus distributed throughout SEA-SC, making it an ideal study group to examine specific biogeographic hypotheses. In addition, the diversity of this genus remains poorly understood, and the phylogenetic relationships among species of Leptolalax and closely related Leptobrachella Smith 1928 remain unclear. Herein, we evaluate species-level diversity based on 48 of the 53 described species from throughout the distribution of Leptolalax. Molecular analyses reveal many undescribed species, mostly in southern China and Indochina. Our well-resolved phylogeny based on multiple nuclear DNA markers shows that Leptolalax is not monophyletic with respect to Leptobrachella and, thus, we assign the former to being a junior synonym of the latter. Similarly, analyses reject monophyly of the two subgenera of Leptolalax. The diversification pattern of the group is complex, involving a high degree of sympatry and prevalence of microendemic species. Northern Sundaland (Borneo) and eastern Indochina (Vietnam) appear to have played pivotal roles as geographical centers of diversification, and paleoclimatic changes and tectonic movements seem to have driven the major divergence of clades. Analyses fail to reject an "upstream" colonization hypothesis, and, thus, the genus appears to have originated in Sundaland and then colonized mainland Asia. Our results reveal that both vicariance and dispersal are responsible for current distribution patterns in the genus.

Mitochondrial genomes uncover the maternal history of the Pamir populations.

The Pamirs, among the world's highest mountains in Central Asia, are one of homelands with the most extreme high altitude for several ethnic groups. The settlement history of modern humans on the Pamirs remains still opaque. Herein, we have sequenced the mitochondrial DNA (mtDNA) genomes of 382 individuals belonging to eight populations from the Pamirs and the surrounding lowlands in Central Asia. We construct the Central Asian (including both highlanders and lowlanders) mtDNA haplogroup tree at the highest resolution. All the matrilineal components are assigned into the defined mtDNA haplogroups in East and West Eurasians. No basal lineages that directly emanate from the Eurasian founder macrohaplogroups M, N, and R are found. Our data support the origin of Central Asian being the result of East-West Eurasian admixture. The coalescence ages for more than 93% mtDNA lineages in Central Asians are dated after the last glacial maximum (LGM). The post-LGM and/or later dispersals/admixtures play dominant roles in shaping the maternal gene pool of Central Asians. More importantly, our analyses reveal the mtDNA heterogeneity in the Pamir highlanders, not only between the Turkic Kyrgyz and the Indo-European Tajik groups, but also among three highland Tajiks. No evidence supports positive selection or relaxation of selective constraints in the mtDNAs of highlanders as compared to that of lowlanders. Our results suggest a complex history for the peopling of Pamirs by multiple waves of migrations from various genetic resources during different time scales.

Out of Southern East Asia of the Brown Rat Revealed by Large-Scale Genome Sequencing.

The geographic origin and migration of the brown rat (Rattus norvegicus) remain subjects of considerable debate. In this study, we sequenced whole genomes of 110 wild brown rats with a diverse world-wide representation. We reveal that brown rats migrated out of southern East Asia, rather than northern Asia as formerly suggested, into the Middle East and then to Europe and Africa, thousands of years ago. Comparison of genomes from different geographical populations reveals that many genes involved in the immune system experienced positive selection in the wild brown rat.

Identity-by-descent refines mapping of candidate regions for preaxial polydactyly II /III in a large Chinese pedigree.

Preaxial polydactyly (PPD) is congenital hand malformation characterized by the duplication of digit. Herein, we scan the genome-wide SNPs for a large Chinese family with PPD-II/III. We employ the refined IBD algorithm to identify the identity-by-decent (IBD) segments and compare the frequency among the patients and normal relatives. A total of 72 markers of 0.01 percentile of the permutation are identified as the peak signals. Among of them, 57markers locate on chromosome 7q36 which is associated with PPD. Further analyses refine the mapping of candidate region in chromosome 7q36 into two 380 Kb fragments within LMBR1 and SHH respectively. IBD approach is a suitable method for mapping causal gene of human disease. Target-enrichment sequencing as well as functional experiments are required to illustrate the pathogenic mechanisms for PPD in the future.

Rapid Evolution of Genes Involved in Learning and Energy Metabolism for Domestication of the Laboratory Rat.

The laboratory rat, widely used in biomedical research, is domesticated from wild brown rat. The origin and genetic mechanism underlying domestication of the laboratory rat remain largely elusive. In the present study, large scale genomes supported a single origin for the laboratory rat, possibly from a sister group to wild rats from Europe/Africa/Middle East. Genomic and transcriptomic analyses uncovered many artificially selected genes (e.g., FOXP2, B3GAT1, and CLOCK) involved in the nervous system. These genes associate with learning ability and regulation of circadian rhythm, which likely enabled the successful domestication of the laboratory rat. Particularly, many genes, including mitochondrial genes responsible for energy metabolism, displayed a substantially increased expression in the brain of laboratory rats compared with wild rats. Our findings demystify the origin and evolution of this model animal, and provide insight into the process of its domestication.

Comparative population genomics reveals genetic basis underlying body size of domestic chickens.

Body size is the most important economic trait for animal production and breeding. Several hundreds of loci have been reported to be associated with growth trait and body weight in chickens. The loci are mapped to large genomic regions due to the low density and limited number of genetic markers in previous studies. Herein, we employed comparative population genomics to identify genetic basis underlying the small body size of Yuanbao chicken (a famous ornamental chicken) based on 89 whole genomes. The most significant signal was mapped to the BMP10 gene, whose expression was upregulated in the Yuanbao chicken. Overexpression of BMP10 induced a significant decrease in body length by inhibiting angiogenic vessel development in zebrafish. In addition, three other loci on chromosomes 1, 2, and 24 were also identified to be potentially involved in the development of body size. Our results provide a paradigm shift in identification of novel loci controlling body size variation, availing a fast and efficient strategy. These loci, particularly BMP10, add insights into ongoing research of the evolution of body size under artificial selection and have important implications for future chicken breeding.

Positive selection rather than relaxation of functional constraint drives the evolution of vision during chicken domestication.

As noted by Darwin, chickens have the greatest phenotypic diversity of all birds, but an interesting evolutionary difference between domestic chickens and their wild ancestor, the Red Junglefowl, is their comparatively weaker vision. Existing theories suggest that diminished visual prowess among domestic chickens reflect changes driven by the relaxation of functional constraints on vision, but the evidence identifying the underlying genetic mechanisms responsible for this change has not been definitively characterized. Here, a genome-wide analysis of the domestic chicken and Red Junglefowl genomes showed significant enrichment for positively selected genes involved in the development of vision. There were significant differences between domestic chickens and their wild ancestors regarding the level of mRNA expression for these genes in the retina. Numerous additional genes involved in the development of vision also showed significant differences in mRNA expression between domestic chickens and their wild ancestors, particularly for genes associated with phototransduction and photoreceptor development, such as RHO (rhodopsin), GUCA1A, PDE6B and NR2E3. Finally, we characterized the potential role of the VIT gene in vision, which experienced positive selection and downregulated expression in the retina of the village chicken. Overall, our results suggest that positive selection, rather than relaxation of purifying selection, contributed to the evolution of vision in domestic chickens. The progenitors of domestic chickens harboring weaker vision may have showed a reduced fear response and vigilance, making them easier to be unconsciously selected and/or domesticated.

An integrated model of clinical information and gene expression for prediction of survival in ovarian cancer patients.

Accumulating evidence shows that clinical factors alone are not adequate for predicting the survival of patients with ovarian cancer (OvCa), and many genes have been found to be associated with OvCa prognosis. The objective of this study was to develop a model that integrates clinical information and a gene signature to predict the survival durations of patients diagnosed with OvCa. We constructed mRNA and microRNA expression profiles and gathered the corresponding clinical data of 552 OvCa patients and 8 normal controls from The Cancer Genome Atlas. Using univariate Cox regression followed by a permutation test, elastic net-regulated Cox regression, and ridge regression, we generated a prognosis index consisting of 2 clinical variables, 7 protective mRNAs, 12 risky mRNAs, and 1 protective microRNA. The area under the curve of the receiver operating characteristic of the integrated clinical-and-gene model was 0.756, larger than that of the clinical-alone model (0.686) or the gene-alone model (0.703). OvCa patients in the high-risk group had a significantly shorter overall survival time compared with patients in the low-risk group (hazard ratio = 8.374, 95% confidence interval = 4.444-15.780, P = 4.90 × 10(-11), by the Wald test). The reliability of the gene signature was confirmed by a public external data set from the Gene Expression Omnibus. Our conclusions that we have identified an integrated clinical-and-gene model superior to the traditional clinical-alone model in ascertaining the survival prognosis of patients with OvCa. Our findings may prove valuable for improving the clinical management of OvCa.

Divergence of dim-light vision among bats (order: Chiroptera) as estimated by molecular and electrophysiological methods.

Dim-light vision is present in all bats, but is divergent among species. Old-World fruit bats (Pteropodidae) have fully developed eyes; the eyes of insectivorous bats are generally degraded, and these bats rely on well-developed echolocation. An exception is the Emballonuridae, which are capable of laryngeal echolocation but prefer to use vision for navigation and have normal eyes. In this study, integrated methods, comprising manganese-enhanced magnetic resonance imaging (MEMRI), f-VEP and RNA-seq, were utilized to verify the divergence. The results of MEMRI showed that Pteropodidae bats have a much larger superior colliculus (SC)/ inferior colliculus (IC) volume ratio (3:1) than insectivorous bats (1:7). Furthermore, the absolute visual thresholds (log cd/m(2)•s) of Pteropodidae (-6.30 and -6.37) and Emballonuridae (-3.71) bats were lower than those of other insectivorous bats (-1.90). Finally, genes related to the visual pathway showed signs of positive selection, convergent evolution, upregulation and similar gene expression patterns in Pteropodidae and Emballonuridae bats. Different results imply that Pteropodidae and Emballonuridae bats have more developed vision than the insectivorous bats and suggest that further research on bat behavior is warranted.